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Models with range-free frustrated Ising spin interaction and additional Hubbard interaction are 
treated exactly by means of the discrete time slicing method of Grassmann field theory. Critical and 
tricritical points, spin- and charge correlations, and the fermion propagator, are derived as a function 
of temperature, chemical potential p, of the Hubbard coupling U, and of the spin glass energy J. U 
is allowed to be either repulsive (U > 0) or attractive (U < 0). Cuts through the multi-dimensional 
phase diagram are obtained. Analytical and numerical evaluations take important replica symmetry 
breaking (_R5'-B)-effects into account. Results for the ordered phase are given at least in one-step 
approximation (1RSB), for T — Owe report the first two-, three-, and four-step calculations (4RSB) 
for fermionic spin glasses. The use of exact relations and invariances under RSB together with 
2i?S'_B-calculations for all fillings and 4i?5'-B-solutions for half filling allow to model exact solutions 
by interpolation. For T — 0, our numerical results provide strong evidence that the exact spin glass 
pseudogap obeys p(E) — ci\E — Ef\ for energies close to the Fermi level with ci ~ 0.13. Rapid 
convergence of p'(Ef) under increasing order of RSB is observed and p"(E) is evaluated to estimate 
subleading powers. Over a wide range of the pseudogap and after a small transient regime p(E) 
regains a linear shape with larger slope and a small S-like perturbation. The leading term resembles 
the Efros-Shklovskii Coulomb pseudogap of two-dimensional localized disordered fermionic systems. 
Beyond half filling we obtain a v — 1 ~ (p — U) 2 ,p > U, dependence of the fermion filling factor 
v. We find a half filling transition between a phase for U > p, where the Fermi level lies inside the 
Hubbard gap, into a phase where p{> U) is located at the center of the upper spin glass pseudogap 
(SG-gap). For p> U the Hubbard gap combines with the lower one of two SG-gaps (phase I), while 
for p < U it joins the sole SG-gap which exists in this half- filling regime (phase II). Shoulders of the 
combined gaps are shaped by RSB due to spin glass order. We predict scaling behaviour at the half 
filling transition which becomes continuous due to ooRSB. Implications of the half-filling transition 
between the deeper insulating phase II and phase I for the eventual derealization by additional 
hopping processes in itinerant model extensions are discussed. Possible metal-insulator transition 
scenarios are described. Generalizations to random Hubbard coupling and alloy models as well as 
frustrated magnetic interactions with ferro- or antiferromagnetic components are also considered 
separately. 



I. INTRODUCTION 



The interplay between complex magnetic order and electronic properties including transport attracts growing inter- 
est despite an already long history of active and successful research in this field. Reviews of many physical phenomena 
which result from competing interactions in frustrated magnetic systems can be found for example in Ref. ^ while 
coupling with electronic transport is discussed in Ref. The particular role of Hubbard interactions in the context 
of glassy order was considered for random transition metal alloys in Ref. || already more than three decades ago. 
This work described the multiple competition between different interactions, which link magnetism with electronic 
transport such as mobile carriers in contact with spin glass order, and noted also the interference of the Kondo effect. 
Moreover an analogy with Anderson localization theory was formulated. At that time the theory of Anderson local- 
ization continued to progress over many years and was finally shaped in its replicated field theoretic version in the 
early eighties, while, independently and almost simultaneously.-Hiith this event, Parisi discovered the highly complex 
replica symmetry broken solution of range-free spin glass orderoa. These theories remained remarkably disconnected 
for a long time, if one disregards for example the analogy discussed in Ref. |^. On one hand, one knew that the 
unitary universality class of Anderson localization takes into account broken time-reversal invariance in the case of 
random magnetic scattering centers. However this class, derived for noninteracting disordered systems, cannot simply 
be assumed to describe localization of charge carriers due to their exposure to spin glass order, which itself is an effect 
of magnetic interaction. While the potential power of random magnetic order to localize charge carriers is evident, 
the specific part played by the many body interactions, which can create gaps in the density of states near the Fermi 
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level and other complicated fluctuation effects, must also be, taken care of. A well-known case where disorder and 
interaction are essential is that of Mott- Anderson transitionso. The role of magnetic fluctuations due to singularities 
in the triplet channel remain a matter of concern. In recent years, thje,Mott transition of the Hubbard model, ig- 
noring the long-range Coulomb interaction, raised a lot of interest tooLrcl. The theory of insulators and in particular 
of antiferromagnetic materials is to a large extent connected with the Hubbard model. The apparent link between 
magnetism and transport properties and the presence of a metal-insulator transition renders the model particularly 
attractive. Far from simple however is the many body theory and even the mean field-theory of the Hubbard models. 
Expansions around its local, so-called atomic, limit have been elaborated by MetzneniJ. The necessity to generate the 
effective magnetic interactions and to calculate their magnetic effects in cooperation with transport properties on the 
basis of the simple atomic limit requires a lot from a perturbation technique. Effective magnetic interaction models 
were also derived, for example the famous t — J model, which project onto a subset of relevant variables. 

An example which combines, on the same footing, the Hubbard interaction with magnetic effects of a (nonlocal) 
spin interaction, was introduced in Ref. 11 and found necessary in order to deal with weak localization in He 3 . 



Our present paper introduces and analyzes a model defined by the grand-canonical Hamiltonian 

U = U ^ nifn-il ~ ^ JijCTjCFj - H^n irJ (1) 

i (ij) ia 

where the (ij) run over all pairs and the random magnetic interaction involves Gaussian distributed independent 
couplings Jij. Spin and charge operators cr^ = — and n = are given in terms of fermion operators c. The 
Hubbard interaction is nonrandom (except in a final section |x|). The model aims at the description of insulating 
phases, where the leading role is played by the cooperation of Hubbard- and frustrated spin-interactions. One can 
view the model as the localized limit of a Hubbard model with frustrated spin interactions and hence, in addition, as 
a starting point with highly nontrivial magnetic behaviour for expansions with respect to electron hopping. We show 
that the model is soluble for range-free Jy to the extent that one can derive the full replica symmetry broken solution. 
The study of such insulating models, or of the insulating side of itinerant models, helps to analyze the precursor of 
the onset of transport-related phenomena which depend on the interplay between magnetism and energy gaps (in 
particular at the Fermi level). The breakdown of magnetism and a subsequent and related filling of the gap, when 
hopping is added, can induce an insulator metal transition. 

Another goal of the research presented in this paper is to study the enhancement of insulating properties and 
the transition between distinguishable types of insulating phases. The enhancement we find here originates in a 
constructive superposition of gap widths. It is in general of considerable technological relevance to identify all 
elementary interaction mechanisms which reinforce or change in a controlled way the insulating tendency of materials 
and, at the same time, to find relationships with magnetic order. Random magnetic interactions deepen insulating 
properties by strong density of states depletion near the Fermi level and in addition by the randomness which generates 
a different type of Anderson localization. 

The present work reports progress in understanding insulating fcrmionic models which incorporate HubbarcLmtpr.- 
actions in the fermionic extension of the range-free Sherrington-Kirkpatrick spin glass interaction (iSlf-model)Bt3E3. 
We note in passing that the model contains as special cases not only the standard SK-mode\ (at imaginary chem- 
ical potential fi = iirT/2) and the atomic limit of the Hubbard model, but also the Ghatak-Sherrington (spin 1) 
modeltJ and of course, for vanishing Hubbard interaction, the fermionic SK-model, which was analyzed in Ref. p"5|on 
li?5B-level. When a strong Hubbard repulsion effectively evacuates double occupied sites the model operates in an 
almost-three-state per site space. In this regime it is also similar to the model of randomly placed strong {/-scattering 
centeralj provided the fermions are localized. 

The article is organized in three main parts A — C and two additional ones D, E. Part A is purely analytical, while 
B and C contain the numerical analyses of the exact selfconsistent equations obtained from A. In particular, the 
reader will find in 

A) [section |]j the derivation of the analytic solution by means of the discrete time slicing method, where all fields 
representing t he Hubb ard interaction are integrated out in perfect coexistence with the spin glass related fields, 



B) [sections III VII |- a detailed numerical analysis of critical and tricritical temperatures, cuts through the multi- 
dimensional phase diagram, solutions of the coupled selfconsistent integral equations in the ordered phase, their 
subsequent use in the calculation of the fermion propagator including the density of states and a calculation of the 
fermion concentration, all given as a function of the chemical potential /i. The full range of /x for which the spin glass 
phase can exist is covered. Particular attention is paid to the large replica symmetry breaking (RSB-) effects at low 
temperatures including T = 0. For T — we present the first 2-step RS ^-evaluation for fermionic spin glasses (with 
arbitrary filli ng) , 



C) [section VIII the 4-step RSB solutions for half filling which lead, by interpolation with known exact relations, to 



a prediction of the exact form of the spin glass pseudogap. The nonanalytic behaviour is determined and a comparison 
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with the Efros-Shklovskii Coulomb pseudogap included. 

In the remaining smaller part D randomness in the Hubbard coupling is included in order to control changes in the 
domain of discontinuous and continuous phase transitions, and finally in E, implications of the gap structure derived 
in part B are discussed. In particular we include in 

D) [section IX] a phenomenological discussion of the onset of diffusive transport beyond a metal-insulator transition. 
Two scenarios for metal-insulator transitions governed by the random and frustrated magnetic interaction of spin 
glasses are described. This part serves to describe one of the major aims of our current research. Finally in 

E) [section X, XI] the effect of several distributions of random Hubbard interaction on the freezing temperature Tf 
are analyzed, and a discussion of ferro- or antiferromagnetic effects induced by nonzero mean values of concludes 
the paper. 

Continued efforts to approach the (best possible) analytical T = 0-solution for the mean field theory of non-itinerant 
spin glasses form a part of the paper. Sufficiently good qualitative, if not quantitative, knowledge of this localized 
T = solution seems necessary, or at least a helpful basis for the analysis of itinerant glassy magnetism in fermionic 
systems. This can easily be understood if one reviews the strong low temperature RS B-effects on relevant properties 
(like shape, effective width, and stability against thermal filling) of gaps in the band structure. Pseudogaps at the 
Fermi level or smallness of the density of states in the vicinity of Ep are features which undoubtedly are relevant for 
two-particle correlation functions such as conductivities. For this reason RSB has the power to influence T — phase 
transitions and must be controlled, in particular with respect to preditions for metal-insulator transitions. It is also 
interesting to compare our results with spin dynamical fluctuation effects in disordered or translationally invariant 
quantum-spin systems. Unfortunately (in the sense of solvability) T = phase transitions often occur as quantum spin 
glass transitions which involve utterly complicated connections between quantum-spin dynamics, critical behaviour 
and RSB. 

Implications of RSB in many body quantum physics have not yet been sufficiently explored. We mention two 
examples: i) localization theory involves non self-averaging quantities and a RSB-free theory is questionable with 
respect to these observables, and ii) recent theories of quantum critical paints in low-dimensional transverse field 
Ising systems have shown the importance of infinite-disorder fixed pointscll, which were expected to be applicable 
to quantum spin glass transitions as well. Their relevance in high dimensions (domain of attraction) and even in 
mean-field limits must be analyzed. Although the comparison between non-replicated and replicated methods is often 
difficult, the claim of relevance of strong coupling theories can well be related to the strong coupling theory needed 
to incorporate complete RSB, which was a motivation of the present research. We cannot answer these questions in 
the present paper, but our attempt to find analytical solutions at T — for the present infinite-range problem is also 
meant as a step towards finite-dimensional systems and the control of BSB-relevance. If one wants to understand the 
term 'relevance' in the renormalization-group sense, it is important that continuous T = phase transitions can exist 
in the system. While the magnetic breakdown is discontinuous in the mean field case and probably in general too, 
we derive and emphasize the special nature of half-filling transitions of the present model. Moreover this transition 
is sharp only at T — and becomes continuous only due to infinite RSB, while in other cases RSB renders a first 
order transition even more discontinuous (one example given in Ref. |l8| ). ■— ■ 

Also recently and during progress of the present work, a numerical approach using infinite RSB equations^ was 
proposed which employs an interesting analogy between the 3SAT-optimization problemE3£3 and the standard SK- 
model at T = 0. In our approach, apart from considering more general model, we follow a different approach which 
attempts to evaluate the maximum feasible number of integrations involved in ooRSB equations. In this paper the 
highest accuracy is so far the 2BSB-approximation of the exact solution for arbitrary values of the chemical potential 
and ARSB for particle hole symmetry /i = U/2 (half filling). 

Since the techniques applied in section [n] are readily extended to include statistical distributions of the Hubbard 
interaction and/or different interaction strength for A and B-type atoms in A-B alloys we extended our analysis to 
study models described by the Hamiltonian 7i + 5TC with 

8H=y j SUjn^nn, (2) 

i 

where 5U obeys a wide variety of different distributions. We evaluate changes in critical temperatures and their 
effect on the phase diagram. This part helps to gain insight into the competition between repulsive and attractive 
interactions on one side and spin glass order on the other. The strong asymmetry between positive and negative U are 
incorporated in one model and weighted by [/-distributions. A distribution-invariant point is found (for symmetric 
distributions), i.e. a point where the critical temperature does not depend on [/-fluctuations at all. 

Due to the absence of spin- and charge-quantum dynamics (in the hopping free model) one may expect the model 
to be more easily solvable than the Hubbard model, the solution of which faces the high obstacle of complicated 
dynamics. However randomness of the magnetic interaction introduces the complication of replica symmetry breaking 
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at all temperatures below freezing and in particular at T = too. The additional introduction of hopping leads to 
further complications, since all problems of the Hubbard model are then contained too. Let us mention the two most 
important classes of itinerant model extensions. These are given by the single- fermion and (fermion)-pair hopping 
Hamiltonians by 

a ) T-Lt = ^2 t ij4* c i<r I b) H p t mr = X]*if c li c iT c iT c j-L' ( 3 ) 

(hi) 

which respectively have the tendency to delocalize single fermions or pairs of fermions. Each of these additional 



hopping Hamiltonians can be combined with models (1) and (2). In section [X we describe qualitatively effects of 
the single-fermion hopping term Tit- A perturbation theory in the fermion hopping elements requires detailed (if not 
exact) knowledge of the fermion propagator in the insulating limit, which is analyzed in section VIIIC. 

The pair- hopping Hamiltonian 7i^ air , on the other hand, can cause a coherent pair state with various consequences. 
Related questions concern the possible existence of microscopic superconducting glasses. An interesting simplification 
of 3b) if compared to 3a) is the fact that H^ air commutes with Hamiltonian ([!]) , which means that the spin glass field 
Q has a static saddle point. Quantum dynamics of the SU(2) generators CfCj,, cjcj does not invalidate a Q-static 
theory (unlike the H t Hamiltonian). Since the pair operators are SU(2) generators as well, a distribution of the 
pair-hopping elements analogous to the Jij in an SK-model would generate an XY superconducting glass model. 
Its description requires a quantum-dynamic order parameter and thus falls into another class of low temperature 
behaviour. r~in 

Experiments on two-dimensional electron gases (2DEG) with strong ferromagnetic coupling to local momentsE!rE3 
demonstrated the relation between transport and magnetism in the presence of disorder and correlation. A spectacular 
giant magnetoresistance was observed in small fields while a transition into a quantum Hall liquid took place at high 
fields. The 2D-geometries were intentionally used to rule out orbital effects and in Ref. |2l] to suppress magnetic 
order. Results of our paper show that a strong ferromagnetic coupling to a spin glass ordered II- VI semiconductor, for 
example, would present a novel case with certain new fundamental questions reaching beyond the mere comparison 
of the difference between paramagnetic- and spin glass response. This experimental situation would test localization 
mechanisms of more than one kind (and different from the unconfirmed^ possibility of polaron formation): given 
that the strong coupling between 2DEG and randomly oriented magnetic moments support Anderson localization 
while a spin glass pseudo-gap, which we discuss in this paper in comparison with the Efros-Shklovskii gap, is expected 
to contribute a particular stretched-exponcntial low temperature resistance. Moreover the spin glass order does not 
smear the correlation-induced level correlation, but increases it further as described below. 

One final remark refers to the energy scale used in our article: since we focus on the spin glass phase generated 
by a Gaussian-distributed random spin interaction the most important energy scale is set by J — \/ ((Jij — (Jij)) 2 }- 
We work hence in J -normalized dimensionless energy variables in order to ease reading. Only when needed, this 
normalization is mentioned explicitly. 



II. EXACT DERIVATION OF SELFCONSISTENT SCHEMES FOR FRUSTRATED RANGE-FREE 
MAGNETIC INTERACTION AND HUBBARD INTERACTION 

Working the double-interaction model (Q) into a replicated Grassmann field theory is a standard procedure. The 
Grassmann field technique is explained for example in Ref. |2^. Usually a continuous time limit is taken before the 
fermionic fields are integrated out. This involves an unknown infinite constant which is removed by a regularization. 
The discrete time-slicing method can instead be used and all integrations performed before the continuum time limit 
is taken at the end. This is known to be the exact procedure of field theory. It avoids singularities of the symbolic 
continuous-time formalism and hence does not require regularization. Its feasibility depends on whether a quadratic 
form in the Grassmann fields can be obtained and in particular be diagonalized. In spite of an 0(ip s ) effective spin 
glass interaction (before decoupling) and the subsequent complication by RSB, we found it convenient to solve the 
problem in this exact way (up to the analysis of RSB which is the only source of approximations in the paper). 

The Hubbard interaction, local in real space and in time, competes with the infinite range SK Hamiltonian and 
its statistical spin correlations, which are local in time but also involve time-independent (infinite-ranged in time) 
magnetic disorder correlations. Both interactions can be decoupled altogether. Not only the Grassmann integrals can 
be performed but also the Hubbard-decoupling fields will be integrated out exactly. This procedure is applied first 



to the partition function in II A and subsequently to the fermion propagator in II B, where it requires a modification 



This amounts to a diagonalization of the problem, exempting only the final solution of infinite-step replica symmetry 



breaking. Reaching a four-step accuracy at T = in section VIII turns out to be sufficient for drawing the physically 
relevant conclusions. 
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A. The disorder averaged grand canonical partition function and thermodynamic behaviour 



The grand canonical partition function of model (jl|) for fermions in a particle bath with chemical potential fi can 



be expressed in terms of the anticommuting eigenvalues ■0", ip a of the fermion operators 



This is achieved by 



inserting the representation of 1 in terms of fermionic coherent states at M — 1 equidistant time slices of the imaginary 
time interval < r < (3 = l/(fcsT). The index a denotes an n-fold replication of the system. Then, as usual in the 
framework of the replica trick, the averaged free energy F = —T(lnZ) — lim n ^o((l — Z n )/n) is calculated from the 
disorder average of 



M-i „ 

Z " = II II / #f CT (r fe )#? CT (r fe ) exp 



exp 



k— i,a,a 

M-2 



M-2 



e E E (^(CfeiJ.c^D+rt-ctroj.ctrM-i))) 



a fe=0 



fe=0 



(4) 



where 



/c a (0, v) = E ^E^^+ i )*^)^'^+ i ) <r '^'(^) 



E 



/<Ec(%i)CK)- 



(5) 



The disorder average, performed over the chosen Gaussian distribution of frustrated magnetic couplings Jjj (but fixed 
Hubbard-coupling U of arbitrary sign) leads to an effective functional with 8-fermion correlation. This integration 
and the subsequent decoupling of the 8-fermion interaction (which is equivalent to the usual 4-spin interaction) are 
standard and not repeated here in detail. We just mention that the Gaussian decoupling by means of matrix fields 
Qi b {Tk,Tk') initially depend on two times as well as replica labels, since those carry along all degrees of freedom as 
obvious from the coupling term Q ab (rk, Tk 1 ) Yla 4'a(' T k)< :T ' l l J a( T k) J2a< ^M T fc') cr VM T fc')- The decoupling fields for the 
Hubbard interaction, expressed above in terms of charge and spin operators, are dynamic; we denote them by ctk and 
7fc for each time instant Tfc. The spin glass Q- fields can in general also show quantum dynamical behaviour. This 
would be the case if transport processes were allowed. As in quantum Heisenberg spin glassesE-3 or e.g. in the famous 
Hubbard model in infinite dimensions^, quantum dynamics poses even a major obstacle in finding an exact analytical 
mean-field solution. We remark that the present techniques in principle allow to include quantum spin- or charge- 
dynamical effects but exact solutions may not be found even in the mean field limit. Here we consider only cases in 
which the mean of is zero or positive so that in the following decomposition 



Qf{Tk,T k ,) = q aa (r k - Tk ,)5 ab + q ab (l - 6 ab ) + SQf(r k , tv). 



(6) 



i b {i~k, Tk>) have a purely static saddle point. Thus one can take advantage of the time 
q, noting in addition that q ab = ((a a )(a b )) dis is always static unless non-equilibrium 



a saddle point matrix without spatial variation applies, and the site i index has been dropped in its elements q. The 
spatial dependence of the-aaddle point solutions must instead be retained in cases like glassy antiferro- or ferrimagnets, 
two-component magnetsEa or similar phases, where regular spatial structures are not removed by averaging over 
disorder. In the absence of quantum spin dynamics, as is the case for the present non-itinerant model, the spin glass 
order parameter fields Qkk 1 = 
independence of q aa (Tk — tv) 
dynamics is taken into account. Moreover contributions from SQ vanish in the thermodynamic limit. It is worth 
noting that also in the finite-range extension of the present model a non-dynamic (Q aa ) is obtained. In contrast to 
Heisenberg spin glasses or metallic spin glasses, where the so-called static approximation of Q aa (t,t') is in principle 
inadequate and perhaps fails to detect important features, the present model imposes a static average of the diagonal 
Q-ficlds - we emphasize that, within the present paper, approximations are only applied in the context of a finite 
number of replica symmetry breaking steps. At present we implicitly keep complete RSB by means of the Parisi 
matrix Q p ar isi and continue to perform exact operations on the replicated partition function, which is obtained after 
averaging (jij) over a Gaussian distribution of the random magnetic interaction as 
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G 

(i>rO 



/G pG pG 
././. 
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(7) 



5 



CM 



The superscript G stands for normalized Gaussian integrals and {r, b r } runs over all spin decoupling fields z r 
on which the effective field i? e // depends. The spin-fields z r r> serve to decouple all (3(?/> 4 )-terms of the form 
(rpa^Q a-fisifyaip) , which are associated with the r-th level block matrices aligned along the diagonal of the Parisi 
Q-matrixO'EI. The crucial tree structure of Parisi- RSB is thus kept track of in the exponent 



oo b r+1 m r+1 /m r 

W^II E r aik [BZ(H eff (y a ,{z^}) 

r—0 b r = 6 r _|_im r +i/m r + (l — m r -)_i/m r ) 



kk' 



(8) 



where YilZ symbolizes a multiple sum over all &o, b\, 62, boo with boundary conditions mo = 1 and moo = n. 
Integration over Grassmann fields tp will readily yield the determinant of B. This matrix is given by 



(BZ) k+1 , k = -(1 + e M + i^Tuj2 a% + a [VTu/2 7 £ + H e ff (y a , {z r })\ ) 
{B™) h k = l and boundary terms B 0jM -i = B M ,M-i 



where 



H eff (y a ,{z r }) = Vq-n Z + EV^+i 



(9) 
(10) 

(11) 



represents the spin glass effective field. The effective field acquires an additional magnetization contribution Jq < a > 
in case one allows a finite mean value of Jo =< Jij >. Truncating the sum over r after the K-th term with q K +i = 
represents the K-step approximation to full RSB, which we call kRSB and 00RSB respectively. The Gaussian integral 
operator representation 

,-G 



dz 



-z 2 /2 



2tt 



(12) 



provides more transparency especially in reading multiple integrals. We use this definition everywhere in the paper. 
For each replica index a, the Grassmann integration leads to the determinant 



detB" 



M-l 

1 + [] [l + e M + {iat + "it) + e aH eff (y a ,{z r }) 



fe=0 



(13) 



which still depends on the z(-spin)-field configuration. Time slicing of the imaginary time interval < r < 8 was 
performed in M — 1 equidistant steps of size e = B/M. The continuum time limit and hence the exact solution will 
result by taking the limit M — > 00 in the end. Let us further abbreviate by use of the shorthand notation 



/G /■(_! I>ij 

detB = nn / / detB 

(1 + [] [l + + Heft) + V^i^l + Ik)} + II I 1 + - Heff) + VeU/2 (*og - 7 g) 
k fc 

II [(! + ^ + *at) 2 - (eH eff + Jdj]2 7 £) 2 ] ) 



(14) 



and extract the contributions which survive in the continuum time limit, expressed by e — > 0,M — > 00, Me = 8. The 
integration over the Hubbard decoupling fields can be executed exactly which yields in the continuum time limit 

,-G 



detB = 2e^ 



coahtffa - j))e-^ 2 + cosh(/3H eff (y a , {z^}) 



(15) 



Using this result together with Eq.^ and Eq.^| in the expression for the free energy per site (respectively the thermo- 
dynamic potential) F = -~T\im n ^ ((Z (n ^) - 1)/ n one finds (after one more integration over the y a -fields) 



IB 1 00 

0F = -x 2 + f (gi - 1)X + j8 2 J2 m r(q 2 r - q 2 r+1 ) ln2 



— lim / ln( 



r Jz r+l 



G i-G 



{ { / { / C({z } m2 /" 11 | m 3/ni2 \m r /m r -i 

J Z2 •* z\ 



C({z r }) = cosh(BH e f f (0, {z r })) + e- f3{u+ x )/2 cosh(B(^ - (7/2)) 



(16) 



(17) 



G 



where \ = @{q ~ <7i) is the so-called single- valley susceptibility. The problem is now solved up to integrations over the 
Gaussian z-spin field distributions; in other words, we have reduced the thermodynamics of the given model to that 
of the standard SK-mode\, which is contained as the special case /x = iirT/2, U = 0, with a different kernel C. The 
field-independent part of C is responsible for the competition between Hubbard coupling, chemical potential, and the 
self-consistent single valley susceptibility, which must be determined self-consistently, is obvious. The i?5-B-effects are 
notoriously difficult to treat. They are of importance for the low T thermodynamic behaviour, and found to be even 
more crucial for the quantum-dynamical fermionic correlations. These are derived from the generating functional 
below. The model is particle hole symmetric (PHS) at fi = U/2, where the spin glass temperature is maximal; 
under PHS Eq.|l6| still shows a U -dependence which is necessary to explain for example the magnetic breakdown for 
sufficiently strong attractive Hubbard interaction (negative U and for arbitrary fi). 



B. Generating functional for quantum-dynamical fermion correlations 

The imaginary time single-fermion propagator (and N-fermion propagators too) 

can be derived by derivatives with respect to the generating fields rj from the generating functional 



$(ry, r)) = In 



(18) 



(19) 



Let us symbolize Gaussian integrals over arbitrary decoupling-fields by Y[ J° j where {X} stands either for {z}-spin 

{ x } 

or {a, 7}-Hubbard fields, and by J T>^ for the Grassmann field integration. Using 



(20) 



we calculated the inverse matrix B 1 in order to obtain the propagator. We choose r > r' with r = Tfc i , T 1 = Tfe, and 
obtain the (disorder-averaged, local) fermion propagator as 
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M->0 



11/ 'ii/ , .,, 
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M-l rjj 
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(21) 



where JV = Yl I z Yl I I a I detB. Splitting the product of the second bracket into one over the interval fej, kf and 

one for the rest, the resulting form allows to solve the Gaussian Hubbard- integrals. Terms of order e 3 / 2 do not survive 
the continuum limit and the exact solution becomes 



r J { z r} a , Jy"' 



_|_ e P(li-<TH e ff) e (n+<rH ef f-U)T jj^f 



(22) 



After Fourier transformation from imaginary times to imaginary frequencies (e„ = (2n + l)7rT) the fermion Green 
function in 1RSB, which illustrates RS B-cEects on the propagator in a relatively simple form, is obtained as 



G r G r G 

ml — 1 



Z2 J Z' 



c(z',z 2 y 



c(z uZ2 y 



e fj(n+aH aff ) + 1 e /3(fi,-aH ef f)f e P(v,-U+aHef f ) + ^ 

ie n + u + aH e f f 



ie„ + fx-U + aH e ff 



(23) 



with normalizing factor C{z\, z 2 ) = 2e l3 ^ +x ^ 2S> C(zi, z 2 ) and C given by Eq.JlT]. The result reduces to the known limits: 
in case of vanishing effective field H e ff = the propagator of the Hubbard model's atomic limit is retrieved] and, for 
U = 0, the fermionic Ising spin glass limit is retainedcJ. The result Eq.^3] displays the [/-correlation induced two- level 
splitting, which is broadened into bands by the statistically distributed spin glass fields H e ff(z). This smearing 
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of delta-peaks into bands might lead to less singular expansions in powers of hopping terms than those expansions 
encountered in the pure Hubbard model. We shall further solve below the y-integration for all T and at T = also 
the z\ -integration in order to get a most explicit analytic form. The remaining numerical task is still not quite easy, 
since all parameters on which the effective field depends and the Parisi block parameter m = m(/x, U,T) must first 
be found as selfconsistent solutions and then used in the calculation of the density of states. As Eq.^3] suggests the 
procedure can be extended to formulate TV-particle propagators. 



III. DEPENDENCE OF FREEZING TEMPERATURE ON REPULSIVE AND ON ATTRACTIVE 

HUBBARD COUPLINGS U 

The temperature Tf at which spin glass freezing occurs depends on the chemical potential fj, and on the Hubbard 
coupling U. The reduction of doubly occupied sites under increasing U > favors magnetic order and thus counteracts 
the effect of a growing deviation from particle hole symmetry, described by the parameter Aphs = |m — \ Thus the 
freezing temperature tends towards the SK- value in the infinite repulsion limit for fixed Ap^s- On the other hand, 
negative (attractive) Hubbard interactions act in the same direction as 'particle pressure' from the particle reservoir: 
both enhance double occupation. At a critical negative U the magnetic phase will break down. 

The derivation of the nonperturbative selfconsistent Tf-equation follows the standard procedure for the fermionic 
generalization of the SK-model. The results of this section identify the regimes of continuous and of discontinuous 
transitions below a tricritical temperature called T/3. 

The selfconsistent equation for the freezing temperature of the infinite range model reads (in units of J) 

T f = l/{l + e- {1+u ^ ( - 2T ^cosh{{^~U/2)/T f )), (24) 
which reduces in the special case U — and fj, = mT/2 to the standard SK-model result. 

The numerical evaluation of Tf for several fixed values of U, including attractive interactions, is displayed in Figure 
0. The Fi gure illustrates the competition between chemical potential and Hubbard repulsion. 

The vertical dashed lines at fj, = U drawn in Fig.|l| will be found later to represent the exact limit of half-filling 
at T = 0, while particle hole symmetry (not shown in the Figure) lies at /1 = \U ■ Figure |l| shows specific cuts 
through the model's multidimensional phase diagram which provide here a complete picture of its dependence on 
chemical potential and Hubbard interaction. Only the first order thermodynamic transition lines Tfi(/i) (defined by 
the crossing of magnetic and nonmagnetic free energies) all resemble the one shown around /1 ~ 10.9 and depend on 
RSB. They are not yet better determined than in 1RSB. The back-swing observed for Tfi — > may well occur at 
the crossover temperature where the next higher RSB-stcp is required. The full i?S'i?-solution is not yet known, but 
is expected to lie closer to the stability limit. Thus Figure [l] determines the phase diagram's /x-dependence apart from 
small RS B-corrections near Tn = 0. In the large [/-limit, the T/-curve approaches the shape of a box with height 1 
and a width given by the value of the chemical potential at the tricritical point. This is obtained in Eq.^ below. 



1 - 
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FIG. 1. Phase boundaries of the solutions which extremize the free energy (|24|), shown as a function of chemical potential 
fi and for a set of fixed Hubbard couplings U = —1.4,-0.5,0.5,3, and 10 (in units of J). Solid bold lines are solutions of 
Eq.(|2^) in the regime of continuous transitions. Dashed bold lines, beneath the dash-dotted bold tricritical line at Tfz = 1/3, 
are solutions of Eq.(p^) too but describe paramagnetic only stability limits preempted by first order transitions. First order 
transition lines, which depend on RSB, axe shown for 1RSB for U — 10. Dashed vertical lines at /i — U show (as a guide for 
the eye) how far /i can exceed the Hubbard coupling U (which is also the half-filling limit at T — 0) before spin glass order 
breaks down. Each curve is symmetric with respect to its particle-hole symmetric point at \x = U/2. 



A. The tricritical temperature T/3 and corresponding chemical potential /j, „ (U) 

The tricritical temperature, included in Figure |l|, appears as the lower limit of second order freezing transitions and 
remains | independent of U. This is readily found by solving Eq.^J together with the condition cosh((^,— U/2)/Tf^) = 
2 exp((l+f) / (2T/3)), where the 0(q 2 )-term of the order parameter equation changes sign. For deviations from particle 
hole symmetry (half-filling) which are larger than a tricritical value |//. 3 — U/2\ given in Eq.(|2^), the order parameter 
can no more vanish continuously. 

A special feature is the breakdown of spin glass order due to sufficiently strong attractive (negative) U couplings. 
As Figure [l] shows, minimum values in the range of negative Hubbard couplings U m i n and of the chemical potential 
exist, below which a spin glass cannot exist. The tricritical chemical potentials are symmetrical with respect to PHS 
at fi = S- and one finds (recalling that fj,,T,U are given in units of J) 

{T f3 ,n f3 (U)-^}={± ±^arccosh(2e^ 1+u ^)} (25) 
For homogeneous attractive Hubbard interaction 

U > U = -1 - \ln2 (26) 
limits the regime of continuous transitions. In a very small range U < U$ discontinuous transitions still exist. 



B. Implications of magnetic breakdown for large negative U < U m in 

The destruction of magnetic order by sufficiently attractive Hubbard interaction is due to a redistribution of density 
of states. A critically high density of paired states breaks down the magnetic order and in principle would create a 
coherent superconducting state in cases when fermion hopping- or pair-hopping mechanisms like those given in Eq.^ 
are present. The nonlocal magnetic interaction J^- by itself can however not generate a coherent superconducting 
state. 

As discussed in a chapter |x|, J-U models with additional random U and in particular with (Ua > 0, Ub < 0) include 
extreme proximity- (vicinity-) effects between local pairs and local moments on neighboring sites. It will be very 
interesting to study the instability of such phases towards derealization by Hamiltonians of the type |^a) and b) . 



IV. BAND- AND GAP-STRUCTURES IN THE PRESENCE OF FRUSTRATED MAGNETIC ORDER 

1. The fermionic density of states of the SK model with Hubbard interaction in replica symmetric approximation 

In models as hard to solve as Hubbard- and spin glass models the density of states has the virtue of being one 
of the simpler quantities with yet a broad physical significance. It is important for many one particle properties. 
Thermodynamic behaviour as observed in quantities like particle number (as a function of chemical potential), internal 
energy, specific heat, thermopower etc. depend on it. Detailed properties of the DOS near the Fermi level influence 
the stability of the insulating phase against fermion hopping processes in itinerant model extensions. Due to the Ward 
identity for charge conservationEj one finds that the DOS also determines the amplitude of the diffusive two-particle 
propagator. The DOS can predict the existence of a metal-insulator transition by means of the filling of an energy 
gap at the Fermi level; however it is insufficient to decide where and under which conditions exactly such a transition 
takes place, since this requires to calculate the decay of the diffusion constant by means of a two-particle correlation 
function. It is one goal of the present paper to construct the basis for such calculations. 

The replica symmetric approximation is known to be unstable everywhere in the ordered phase (in fermionic extensions 







as well as in standard SK- models). While i^SB-corrections remain small only down to characteristic crossover 
temperatures, some crude features like the splitting of bands remain correct even down to zero temperature. It 
provides thus a simple and partially meaningful basis for improved calculations. In particular, under its simplifications 
in the zero temperature limit, the analytical expression for the Green function's spectral weight p(e) = — -^Im{G R (e)} 
is simple and reads in terms of energy variable E = e + p as 

p (E) ee p(E)\ 0RSB = -±= e ^ E+ ^ 2 /^d( mi n{-x, U-p- x/2] - E) (27) 

+ e -(E-Uf/(2q) ( E _(2U-(x + x/2))8(n - x/2 -E) + e^"^"*/ 2 ) 2 /^)^ - max[U + X , M + x/2]), 

where the parameters must be understood as q = q(T = 0, p, U) \ 0RSB , x = x(T = 0, p, U)\ 0RSB = %(T = 0, p, U) \ 0RSB ■ 
These parameters q and x are determined by the zero temperature limit of their selfconsistency equations given by 

q = 1 - \v ~ 1|, u-l = erf( ^ U _X /2 )e(p - § - U), X = ,[T e (,-u-x/2f/(2 g)9{fl - /2 _ u} (28) 

y/2q 2 \J nq 



Eq.^ correctly predicts the crude band structure: 

// U < ^ - x/2 :-»■ three bands for E < U — p — x/2, 2U - p + x/2 < E < p - x/2, E > p + x/2 
If U > p — x/2 :— > two bands for E < — x, E > U + x 

The values of these band limits (at T = 0) will change under improvement of the approximation. That means, with 
each step of replica symmetry breaking the single valley susceptibility x w ih decrease further towards zero and thus 
deviate more and more from the equilibrium susceptibility. This turns the band limits of the lower approximations 
into crossover lines of the higher ones. Finally these crossover lines will lie dense giving way to a totally softened decay 
of the density of states towards its pseudogap (or Hubbard gap, respectively). In order to avoid misunderstandings we 
remark that the Gaussian bands given in Eq.(|28|) will not shrink into delta peaks when the order parameter vanishes 
at and above Tf. The random interaction maintains the Gaussian broadening above Tf. 

The full Green function (in the finite temperature technique) is given by the spectral representation Q{e{) — 
de ^'^J where e; = (21 + l)nT and the finite T replica symmetric DOS p (E,T) is given by 



■ly/T r°° 



^VTx J-oo 

+ ( e (M-tO/T + e {E-U)/ T) ^ dze -^ 2 e ~(V«z+E-Uf/(2Tx) r l {p _ ^ ^ g)) (2Q) 



with 



CO*, U, T, z) := e-VWMcoahiZ-jyi) + cosh(^-) (30) 



2. Fermion concentration 



The fermion concentration for given chemical potential is calculated using Eq.(|30|) in 

KM, U, T) = £^ dE e{E _^ )/T + i p(E) (31) 

The i?S'-approximation of Eqn.^l] already predicts correctly the existence of a continuous half-filling transition. The 
half-filling transition is defined as a transition from filling v = 1 to v ^ 1. Figure || displays its thermal smearing 
for p = 0.7. The critical point at U c = p — x(p,U c )/2 will be strongly corrected by i?5B-effects (as derived in 
the following sections). Since the single- valley susceptibility decreases with RSB and must be expected to vanish in 
ooRSB the half-filling transition occurs at U c = p(= 0.7 in this case). 
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FIG. 2. Thermal smearing of the sharp half filling transition at T = as shown by the fermion concentration vq(T) for fi = 0.7 
in replica symmetric approximation. Temperatures are given in units of J. The critical Hubbard coupling U c = fj, — U c )/2 
in O-RSB-approximation and its exact position (dot) at fj, = U for vanishing nonequilibrium susceptibility \ (see below). 



One aspect of RSB and of the unbroken approximation becomes obvious from Fig.y: temperatures above 0.15[J] 
are enough to render i?55-cffccts small enough and results obtained in ORSB free from obvious failure; the effects 
become large near T = and better approximations are required. Yet the O-RS'-B-approximation has the virtue to 
predict correctly the existence of a transition and to provide crude analytical solutions (without too many integrals) . 
Moreover, RSB has the power to generate nontrivial critical behaviour in the range-free model; mean field critical 
behaviour only refers to the e limination of critical fluctuations in real space. We shall see this further below from the 
4i?5B-results of section IVIII. 



We conclude this section with a demonstration of the typical difference between band asymmetries generated by 
ferromagnetic order (assuming here a model with nonvanishing ferromagnetic mean interaction (J) > and U = 0, 
Figure^ (left)) and the effects of the Hubbard interaction with U > while (J) = in FigjH (right). We chose a strong 
doping case with a pronounced central band. The spin glass gaps are present in both cases. Ferromagnetic order 
shifts spectral weight from the upper magnetic band into the lower one and simultaneously renders the nonmagnetic 
central band asymmetric as well. On the right side the asymmetry above freezing is only due to the deviation from 
half-filling, which below Tf becomes accentuated by the spin glass gap at fi, which splits off the central band from 
the upper magnetic band. 



Figure 3 appended (fig3.jpg) at e.o.f. 



FIG. 3. Comparison of strong band asymmetry effects induced either by partially ferromagnetic interaction (left) with 
(Jij) = 1.5, n — 0.7, U = and selfconsistent ferromagnetic order or (right Figure) due to a large chemical potential /i = 5 with 
Hubbard coupling U = 4.4 and (Jij) = above half filling. In the latter case the lower spin glass gap is immersed in a large 
Hubbard gap and the upper spin glass gap hosts the Fermi level, while in the left Figure only spin glass gaps, symmetrically 
placed with respect to E — 0, are excavated in the presence of strong additional ferromagnetic order. The high temperature 
regime of the right hand side Figure shows two energy levels with Hubbard splitting in the clean atomic limit; the random 
magnetic interaction with variance J (set equal to 1) distributes their weight into bands which remain well separated for the 
chosen ratio U = 4.4. 



V. SOLUTIONS WITH 1-STEP REPLICA SYMMETRY BREAKING 



Since i?S , B-effects are strongest in the low temperature regime, we show first the 1RSB results of all relevant phys- 
ical parameters of our present J — U model (1) in the range between T = and moderate temperatures T m 0.25( J). 
These are collected in Fig.||. For T = the self-consistency equations show that the (T = 0, U = 0)-solution can be 
used to find the (T = 0, U > 0) solutions simply by applying the shift \i — * fi + U. This relation can be concluded 
from the exponential divergence ~ 6 -/ 3 (m-c / +x/2) ( m case f^g^ ^ > (7/2) of expression (17) as {3 — » oo and, of course, 
is confirmed by our numerical evaluations. 

Figure 4 (fig4.jpg) appended at e.o.f. 
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FIG. 4. liJSB-solutions are shown for fixed Hubbard-coupling U = 0.5J in the chemical potential range from particle-hole 
symmetry, /j, = (7/2, to almost the discontinuous breakdown of the magnetic phase at fj, ~ 1.4, and in the temperature-interval 
< T < j: the T-normalized Parisi parameter a(p,T) = m(p,T)/T, the order parameters qi(p,T), q2(p,T) and their 
difference, the equilibrium susceptibility \{Pi T) and single valley susceptibility Traces of the nonanalytic behaviour 

caused by the half-filling transition at (/i = U + x/2, T — 0) disappear with increasing temperature. 

Using these solutions we obtain the density of states. Characteristic cases are displayed in Fig.^. On the left p{E, T) 
for p = 1, U = 0.5 demonstrates for large filling how the broader Hubbard-SG gap and the upper SG-gap are filled 
by thermal excitations. The right Figure p{E, p) shows a clear cut nonmagnetic central band together with a visible 
combination of gap widths at the very low temperature T = 0.01 and U —\. 

Figure 5 (fig5.jpg) appended at e.o.f. 

FIG. 5. Density of states p(E = e+p, T) for p = 1, U = § (left Figure) showing thermal smearing of Hubbard- and spin glass 
gaps and (in the right Figure) weakly thermalized gaps in p(E,p) for T = 0.01, U = | (all energies in units of J). The position 
of the Fermi level is indicated by a thick line, showing the crossover from a spin glass gap into the combined Hubbard-and spin 
glass gap. Both figures are obtained in 17?S'_B-approximation. 

The symmetry between lower and upper magnetic band demonstrate that the half-filling transition coincides with 
the disappearance of the central band as the chemical potential falls below the critical value. As a guide to the eye 
a thick line shows the position of the Fermi level, which leaves the SG-gap below a critical p(U) and moves deeper 
into the Hubbard regime as U/p further increases. The presence of a central nonmagnetic band is synonymous with 
a deviation from half filling in this model. It exists for chemical potentials within the interval p^ < p < p c , where 
Pk — > for k — oo steps of RSB and p c denotes the threshold value for the spin glass phase. 

VI. HALF FILLING PHASE TRANSITION AT T = 0: EFFECT OF THE COMPETITION BETWEEN 
HUBBARD REPULSION AND CHEMICAL POTENTIAL 

A. Splitting and recombination of gaps 

The half- filling transition is sharply defined at T = 0. We summarize the superposition and splitting of Hubbard- 
and spin glass gap(s) in the following Figure [| The spin glass gap (finite in 0RSB) turns into a pseudogap at the 
Fermi level in the exact solution, since RSB fills spectral weight into the bright/vacant regions exempting only the 
Hubbard gap (blue) and the pseudogap line along E — Ep. Thus Figure |fj exhibits a change between two very 
different insulating phases: in the first phase, for p > U, the Hubbard interaction has effects only remote from the 
Fermi level, while in the half-filled regime (p < U) the Fermi level changes over into the Hubbard gap regime. 



12 




2 



E 



-0.5 



-1.5 



0.2 



0.4 



0.6 



0.8 



i u 



FIG. 6. This Figure illustrates the crossover (along the dashed line) of the Fermi level between different gaps in the 
E, [/-plane at T = 0. Dark grey regions show nonvanishing density of states in O-RSB-approximation, light grey the extension 
of these domains in 1RSB. The DOS-landscape is viewed from above. For a chosen Ef = 0.7 in the example shown), Ef 
pins the upper spin glass gap for U < fj. and moves into the Hubbard gap for U > [i. At T = the central nonmagnetic band 
becomes fully suppressed for any Hubbard repulsion stronger than U = p — x/2, which then implies half filling. In the exact 
solution the gap region is reduced and consists of Hubbard gap and spin glass pseudogap (of zero width). 

The way in which 2RSB and higher corrections fill vacant regions of Figure [6] with small but finite spectral weight 
becomes obvious from the results of our paper. Here we want to avoid information overflow in the Figure. For the 
same reason, small deviations from straight lines are not shown, which occur beyond half filling (but, as found below 
in the paper) are identified as RS B-artifacts. The exact solution finally realizes p(E) — only along the straight line 
at Ep and in the Hubbard gap region. The transition between two different types of insulating phases remains as 
described by Figure ||. 

Another type of soft filling into the regions, where Figure ^| still shows evacuated spectral weight, occurs due to 
thermal effects. The transition between Hubbard/spin glass gap phase into the pure spin glass gap phase will be 
manifested for example in the specific heat (as usual specific heat results are reasonable down to temperatures of the 
order T > 7\. E g k, where E g k denotes the gap energy in kRSB). 

In which way transport processes will scatter weight into the gap regimes is a matter of speculation for the moment. 
Essential parts of the present paper are however concerned with a sufficiently good determination of the fermion 
propagators for the insulating phase, which then provide the means to analyze transport effects in a perturbation 
theory. 

VII. DEPENDENCE OF ORDER PARAMETER, SUSCEPTIBILITIES, AND FILLING FACTOR ON 
THE CHEMICAL POTENTIAL: A 27?S'B-ANALYSIS AT ZERO TEMPERATURE 

The half filling transition requires a more refined analysis of replica symmetry breaking effects. The dependence on 
the chemical potential is therefore calculated in a selfconsistent 2i?S*i?-approximation for the full range of p for which 
the spin glass phase exists. We determine all relevant parameters in this section, distinguish i?5-B-artifacts from true 
/^-dependence, and in particular determine the v(p) which helps to change arbitrarily between dependence on filling 
(fermion concentration) or on chemical potential. The filling factor is also linked to the integrated density of states 
(provides internal consistency check), which is determined in greater detail later in the paper. 
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A. Analytical input for the numerical solution at T = 



The extremization of the free energy at finite temperatures is too complicated to have yet led to an analytical 
solution. The T = 0-limit however allows one to reduce the number of integrations. This facilitates a low T-expansion 
and eases the numerical evaluation. We take advantage of this simplification and derive the complete set of T = 
quantities. The T — ► 0-limit leads to the following energy F(0) (per site and in units of J) which, unlike the thermal 
F(T), depends on three (instead of four) order parameters gi.2,3, two Parisi parameters 02.3 = rri2,a/T, and the 
so-called single valley susceptibility \. Extremizing 



1 _ , _ \ 1 ,9 9n 1 ,9 9s 1 



F(T = 0) = -x(qi - 1) + -Mil - & + ~Ml2 ~ 4) - - / In C (z 2 , z 3 ) (32) 



G /.G 



^3 



where 



C (z 2 ,z 3 ) = i V c aifA g.//(Q.o.».*0+°i(«i-p)/2) [1 er/( — 6 ^ ~ A ge//( °' °' Z2 ' Zs) ~ " l(gl ~ g2) ] 
2 a=±i V 2 (<H - 92) 

+ l efc *„ 8(44 y A r _ /( Ag/ie(^)-g e //(0,0,^,Z3) )j 5fi = fj, — U — x/2 (33) 

2 V v %i - 12) 

with H e ff(0, 0, Z2, 23) = V92 — 93 ^2 + y/ga Z3 leads to six coupled (double-)integral equations 

{d ai , d a2 , d qi , d Q2 , 9 ?3 , ajF(o) - (34) 

where d x means J| (Eq|3^ holds for the (j, — U/2> 0-side particle hole symmetry and up to irrelevant constants). We 
determine numerically all parameters for <5/i > from the T = selfconsistency equations which result from Eq|34[ 
The results for [7 = and half-filling are shown in the Figures below and all finite-?/ results can be obtained from 
it by the transformation ji — ► fi + U for positive /i. We shall use the 2i?5S-solutions to sort out artifacts of the 
approximation and understand the true underlying physics. 

Figs.^ and || display new singular features in 2RSB which, by the help of a comparison with the 0- and li?,!5i?-results, 
can be suspected to be artifacts of the approximations; they nevertheless provide an improved insight into the shape 
of the exact Parisi solution in 00RSB. The analytic T = 0-equations explain the origin of k wiggles for fc-stepped 
RSB. The surprising appearance of singular wiggles at smaller ~ § involves multi-valued solutions and, similarly 
to what happens in first order phase transitions, would imply discontinuous behaviour. Finally one arrives at the 
conclusion that the wiggly behaviour must be associated with - and in fact is caused by - the stepped approximation 
of the continuous Parisi function in replica space. A physical meaning of these solutions could yet remain due to the 
fact that the 00RSB equilibrium solution is approached very slowly in time: finite-step RS B-solutions marked by the 
single- and two-valley susceptibilities might be relevant for non-equilibrium behaviour. 
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FIG. 7. On the left side the Figure shows the selfconsistent solutions for T-normalized Parisi parameters a 1 ,a2 and for the 
order parameters qi, 52, <?3 (5 = 1 not shown) at 2i?5'-B-level as a function of the chemical potential /1; on the right the fermion 
filling factor u(fi) is displayed, where labeled arrow-heads indicate the order of the approximation and the gap-size related 
limits of the half-filling regime at \n\ < x (k= °' 1 ' 2) (f J ')/' 2 - The sma11 

region of a multi- valued filling factor is zoomed in the inset. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

FIG. 8. The equilibrium susceptibility \ an d single- valley susceptibility \ are shown in ORSB, 1RSB, and 2RSB approxima- 
tion at zero temperature. The constant parts extend over the half-filling regime (at T — 0) up to fj, hf given by x/2. 



This discontinuous behaviour will be removed by thermal smearing in 2RSB down to much lower temperatures 
than observed in Fig.|| or in the IRSB-g&p. The latter was seen to produce a smooth function like the expected exact 
solution. Throughout the paper we have chosen the chemical potential to be the independent variable. In order to 
see whether the /x-dependence is particularly susceptible to the soft steps, we inverted v(fJ.) into n{v) and obtained 
the physical quantities as a function of the fermion concentration v as independent variable. It turns out that the 
wiggles seen in v(fi) itself do not remove (compensate) these features from the order parameters and susceptibilities 
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FIG. 9. Order parameters and susceptibilities in as a function of fermion concentration v in 27?S f B-approximation. Vertical 
dotted lines locate wiggly behaviour generated by the Parisi steps in replica space. 

Although some unphysical features are still contained in the T = 0-solutions at any finite RSB, 2i?S*i?-solutions 
are much closer to the exact ones. 

Since the Parisi function q{x) is constructed by piecewise continuous/constant functions of x (0- functions below a 
break-point), the present results indicate that this dependence in replica space transforms either into strong variation 
or into jumps which resemble first order phase transitions. Their positions on the /i-axis scale with the step positions 
along the replica axis (replica dimension). A real discontinuous phase transition occurs however at the magnetic 
breakdown value fi w 0.88 for U — and shifted by U as shown in Figure ^J. The critical filling is independent of the 
Hubbard coupling and roughly given by 1.45. 

By inspection of the selfconsistent equations one expects that the number of wiggles will increase with the number k 
of RSB-steps; the size of the jumps will diminish, thus improve the quality of the approximation, and finally disappear 
in the k — > oo-limit. 



B. The order parameter as a function of chemical potential and T = 0-measure of replica space 

The 2i?S'-B-results for the discretized order parameter and the intervals aj, of their step heights given in Figure 
can be combined and represented in one Figure. An instructive representation of the generalized Parisi-function then 
results in the space of /i and of a variable called a x . The latter replaces the variable x used by Parisi to represent the 
solution q(x). In Parisi's papers, the continuous parameter x was introduced to host all matrix block size parameters 
TOfc. These are found within an interval < x < 1 in the replica limit (in addition they lie dense in the ooRSB 
limit). The function q(x) is discretized by means of ^-function (steps) in any finite order of RSB. In case of the low 
temperature limit, the selfconsistent solutions for matrix block size parameter mt turned out to vanish linearly with 
temperature, which required the use of a rescaled parameter afc = rrifc/T. These parameter values are now hosted in 
an unbounded interval and we consider them as discrete values on an a^-axis, where a x replaces the low temperature 
limit of x/T in a T = calculation. Without temperature-rescaling the rr-interval below the break point, where the 
nontrivial part of the Parisi function shows up, would collapses to a point with q(0) infinitely-valued. The nontrivial 
part of the Parisi function still exists at T = and is well represented on the a^-axis. As a function of a x and of the 
chemical potential we obtain the 2i?5'i?-result of q(a x , fi, T — 0) as given by Fig.|To| 

Figures 10 and 11 appended (figlO-ll.jpg) at e.o.f. 

FIG. 10. 27?S f B-approximation of the T — Parisi order parameter as a function of a x = x/T\t=o in analogy with 
a = lim T ^ m(T) /T and fi. The strongly /i-dependent interval lengths over which the order parameter assumes one of its 
possible three constant values are shown at the /j, = 0-side of the plot. At the large fi ~ 0.88-side the first order regime of 
magnetic breakdown is shown. 
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FIG. 11. The calculated T = limit of the free energy in 2RSB. For < /i < x/2 the half-filling energy is higher (dashed 
line). The inset show the multi- valued regime indicated by an arrow in the large-scale energy. 

At T = it seems possible that the solutions as a function of fi for very high k (close to the exact one at k = oo) 
will be made up of piecewise continuous functions of /i, where discontinuities at the 0(k) wiggles and the intervals 
of continuous behaviour both approach zero. Some resemblance with level crossing/repulsion appears. The wiggles 
must finally be seen to mark a soft stepping of the order parameter in contrast to (and as a consequence of) the hard 
steps along a which are imposed by the definition of the Parisi scheme. One may be tempted to ask for a generalized 
scheme that softens the steps also in the replica dimension before the exact limit is reached. This seems even more 
desirable since the higher order ^SB-calculations below show an extreme flatness of the energy landscape. 
Using the above solutions we calculate the T = 0-limit of the free energy of Eq.(|3^). The solution shown in Fig.[ll| 
confirms the half-filling limit at ji = x/2 and a discontinuous transition near [i « 0.33. Thus it is clear that the 
step-approximation of the order parameter function in replica space leads either to strong variation or to jumps, 
which appear in the form of artificial first order transitions at certain values of the chemical potential; these locations 
correspond to and scale with the positions of the Parisi-steps. No such transition can be expected in the exact solution. 
One could speculate whether [i = is the accumulation point of these discontinuities and thus may be the source of 
special critical behaviour in the [i = point in ooRSB. In this case infinitely many wiggles could render the Parisi 
function a mathematically delicate object at half-filling. While Parisi's finite temperature order parameter function 
q(x) is defined on an interval < x < 1 the present replica interval for T — is unbounded and there is not yet an 
indication that the discontinuities must lie dense in the limit of infinitely many i?S*i?-steps. 

The conclusion of this chapter is that the ^-dependence seen in the susceptibilities will disappear and yield \ — 1, x — 
0, while the order parameter function decreases smoothly asu approaches the magnetic breakdown value and steps 
along the et^-axis must be continuized. In this sense Figur e [10| offers a crude picture which calls for an improved and 



refined treatment. Before this is achieved in section VIII the fermion filling v(n), which is related for T = to the 



largest order parameter qi and also to the integrated density of states, is obtained. 



C. The ^(/independence from 2i?£73-calculation for the whole spin glass regime 

The spin glass regime exists at T = in the range < |/x — U/2\ < fi c « 0.9 + U/2 for positive U. Comparing the 
/ci?S'£?-results for k = 0, 1, 2, 3, 4 together with the known exact value at /J, = 0, T = the regime of validity and the 
quality of 2RSB can be estimated. An interpolation helps to conjecture the form of the expected exact solution^(/x). 
The numerical result can also be used to generate an analytical fit, which in the best case would be the exact solution. 
If not exact, it can still represent an analytical form for the propagator needed in diagram expansions of itinerant 



models for example. Figure 12 shows such an interpolation. The fermion concentration v = dep(e) at T = 
equals 2 - q x = 2 - q aa . 
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Only the zoomed inset in Figure [lj shows significant deviations between the z/(/x) fit-function and the 2RSB- 
calculated points. Viewed more globally over the whole /^-interval a simple quadratic behaviour with small quartic 
corrections fits nicely. The exact behaviour of v(n) appears to be simple. The fit function, which interpolates between 
Vexactin = 0) = 1 and 2RSB- values for \i > 0.12, is close to 1 + ^/i 2 + 0(/i 4 ). The way in which the finite half 
filling limits (y = 1) move to smaller chemical potentials is determined in the following section VIII, confirming the 
interpretation v ~ /.t 2 . 



VIII. THE SHAPE OF THE SPIN GLASS PSEUDO-GAP, ITS COMBINATION WITH THE HUBBARD 
GAP, AND A COMPARISON WITH THE EFROS-SHKLOVSKII COULOMB PSEUDOGAP: 
4RSB-CALCULATION AT HALF FILLING, EXACT RELATIONS 

We extend the above-mentioned extremization of the free energy, Eq.[34|, to a 4i?S*i?-solution. On this level and for 
arbitrary filling, the set of ten coupled integral equations, each having seven coupled integrations, 

{d ai , d a2 , d a3 , d a4 , d qi , d q2 , d q3 , d qi , d q5 , d^}F(0) = (35) 

must be solved selfconsistently. We performed this calculation for half filling where the numerical efforts are reduced 
to the solution of eight coupled equations, since qi = 1 and the x-equation decouples from the rest. If the points in 
replica space, where the order parameter values are sought, would not be required to minimize the free energy too one 
could include several higher RS i?-orders with similar effort. This seems however unnecessary for the present purpose 
and hence was discarded. 



A. 4-step RSB approximation for T = 0-order parameter function and -susceptibilities at half filling 



While a Parisi order parameter function for finite temperatures is defined on an interval < x < 1, the T = order 
parameter function is represented by the temperature rescaled quantity a x = lim^^o x/T, which is the natural choice 
due to the m(T) ~ T dependence observed in our selfconsistcnt calculations. An upper bound for the selfconsistent 
a x -solutions does not seem to exist. The maximum value in ARSB is 5.5 and grows rapidly with the order of RSB. A 
mapping a x — > a x = a x /(l + a x ) has the virtue of re-converting even the T = interval to < a x < 1. The following 
Figures show the obtained selfconsistent values and the fitting continuous order parameter function at T = 0. The 
fitting analytical form indicates that a distribution of error functions is probably involved in the exact analytical 
solution of the problem. 
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FIG. 13. A T = order parameter (analytic model) function, which fits well the calculated 5-step (iRSB) approximation, is 
shown. The upper Figure displays {q — q\ — 1, q2, <?s} and {ai, 0,4}, all obtained selfconsistently, for the original variable 
a x . Below, the same function is shown as a function of a x after mapping to the interval [0..1] by means of a x — > a x /(l + a x ). 



A break point, which is a specific in finite temperature Parisi order parameter function q(x), appears to be absent 
in accordance with the expectation: the break point scales like Xb. P . ~ VT and since a x ~ x/T one expects this point 
to move to infinity in the zero temperature limit. 
The equilibrium and nonequilibrium susceptibilities, obtained from our 4i?S'-B-calculation are shown in Figure 14 
The x-decay is well modeled by a function which decays like 1/k when the number of k RS B-steps goes to infinity. 
Logarithmic corrections or a slower decay cannot yet be ruled out. 




RSB-order 




RSB-order 



FIG. 14. The decay of the single-valley (nonequilibrium) susceptibility \ is shown as a function of interpolated RSB-order; 
calculated points are shown by dots and extrapolation to higher order is possible from a simple and good analytical model fit. 
The growth of the equilibrium susceptibility \ towards 1 is shown in comparison with the x~fh m the Figure below. 



B. 4i?S'B-approximation of the spin glass pseudogap 

Thermal effects cover unphysical remainders of replica symmetry and, as one increases the order of RSB, the 
reappearance of incorrect features shifts to ever smaller temperatures. In a comparable way one needs higher orders 
of RSB to describe exactly smaller and smaller excitation energies above (and below) the Fermi level. In other words, 
e = \E — Ep\ scales with T and the larger the excitation energy is the smaller is the effect of the 'fine structure' of 
the Parisi function for small a x as shown in Figure [l3]. The effect can be seen in the fermionic density of states p at 
T = 0. It will also enter the quantum-dynamical fermion propagator, which is determined by its spectral weight p. 
We improved earlier one-step i?5*i?-calculationst3 by three orders and find that four-step RSB (4RSB) is sufficient 
to predict the shape of the pseudogap. 

We derive the exact 4i?S , B-formula for the fermionic density of states p(E = e + p) for T = as 
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p{E) = e 



y / 2n(q 1 - q 2 ) J 5 



AA 2 (a 4 /a 3 ) [m 2 (1 -04/03) 03/03) / e -( B -x-*.//)V(2(9i-«.))<7« a /oi-i 

J4 J2 

(36) 



using the integral operator J" 6 = = dzb exp(— z^/2)/\/2n acting on the accumulated normalizing factors 



{M(x) 





— X 








0.3/0,2 


(J a 2 /0>l 




, M 2 [x) = 


[X 


Jl 1 





(37) 



and C as given by Eq.|33| after replacement of the 2RSB by the 4i?5B-effective field H e f j . Our numerical evaluation 
of Eqj3q for half-filling and U — is shown in Fig.fL5[ 




FIG. 15. The magnetic band structure at T = and for half- filling is shown for comparison in the replica-symmetric (ORSB), 
1-, 2-, 3-, and 4-step replica symmetry breaking approximation (solid curves labeled with 0, 1, 2, 3, 4 correspondingly). The gap 
energies decrease with the increasing RSB-order and the shape of the gap emerges. The wide scale upper Figure shows the 
decrease of width and position of energy regimes, where the next higher order of RSB bifurcates. 

The bottom Figure zooms the gap regime near Ef- The dashed gap-edge-line 'GEL' shows the exact lines connecting p(E g ) (E g 
denoting the gap edge energies) while the asymptotic power law in the limit e = is indicated by the dashed line p'(e = 0)|e| with 
p'(0) approximated by its 4_RS'_B-value 0.13. An interpolation IPMdos of the density of states between the tiny gap-bottom 
regime following this law and the wide gap regime (where 4RSB is almost exact) with a slope closer to 0.3 is modeled by 
extrapolating a p"-calculation to 00RSB. 
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The dash-dotted straight lines locate the shoulder- height p(E g ^) at the gap edges E Qtk of fc-step RSB: the ratio 
between these heights and the gap-width is invariant (= 1) under a change of the order k of RSB, ie i^S-B-invariant. 
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FIG. 16. First derivative of the density of states p(e), e = E — Ef from replica-symmetric to four-step RSB approximation. 
Gap edges are indicated by points for each approximation and an extrapolation beyond these points towards the Fermi level 
(e = 0) is shown. These lines demonstrate the rapid convergence: lines from 2— ,3—, and 4RSB lines meet in the point 
(0,0.13 + O(lCT 3 )). 
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FIG. 17. The second derivatives of the density of states does not show rapid convergence up to fourth order replica symmetry 
breaking. Dots mark again the gap edge positions e = E g . k in k-th order RSB. Their position is well fitted to the shown curve 
(dashed) p"(E Bth ) = 2.1 - + 7.6E g , k . This serves as a crude estimation of p"(0) in ooRSB. A divergence of p" can 

however not be excluded, which would imply the subleading exponent of p(e) to be smaller than 2. 
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In Figure [h] the gap edge points E g k must reach e = in the ooRSB limit; whether these values E g k for further 
increasing values of k eventually fall on a curve with divergent slope for e — > can not yet be concluded from Figure 
[L7L This affects only the subleading behaviour which determines the curvature within the very small asymptotic 



scaling regime. Figure 18 illustrates the remaining small corrections beyond 4RSB. The interpolating approximation 
to the exact solution shows the smallness of the scaling regime, while the larger linear regime up to e w 1 can be 
viewed as subcritical. Without a closer look it is hardly possible to see that the extrapolation of the almost linear 
DOS-decay would not end in e = but leave a tiny finite fake-gap. This is prevented by the crossover to the scaling 
regime as described in the zoomed part of Figure |l5| . 

Figure 18 (figl8.gif) appended at e.o.f. 



FIG. 18. The magnetic band structure at T = 0, for U — 0, and half-filling. The 4J?S'_B-approximation, which contains 
only a tiny gap, is shown in comparison with the interpolated model (IPMdos) of full RSB matching p(Ef) = and the 
p'(Ef) |e|-law, concluded from RSB-convergence towards p'(Ef) ~ 0.13 of Figure with the regime |e| > 0.3 where the 
"l^RSS-values are effectively exact. 



Effects of the Hubbard interaction: At zero temperature and for half-filling the effect of the Hubbard coupling is to 
spread the spin glass gap by the U and to shift the symmetry point \i — > /i — {7/2. The repulsive Hubbard interaction 
preserves half-filling and maintains the gap even in the limit of ooRSB. 



Figure 19 (figl9.gif) appended at e.o.f. 



FIG. 19. The density of states for vanishing Hubbard interaction and asymptotic e|-behaviour is shown in comparison with 
one for finite Hubbard repulsion U/J = 1. At T = and for half filling the Hubbard gap splits the spin glass pseudogap in two 
halfs and preserves the form of the linear p-decay relative to the gap. 

It is perhaps unexpected to evoke universal critical behaviour in the context of a pseudogap in a range-free (or 
infinite-dimensional) model. However the results indicate that critical correlation equivalents appear to exist in replica 
space. It is in this space that a theory perhaps of renormalization group character should be developed in order to 
analyze the critical behaviour (shape) of the pseudogap. In the better known case of finite range interactions (in finite 
dimensions) the renormalization group approach as defined and applied to the Coulomb pseudo-gap by Johnson and 
KhmclnitskiiEil is a related very interesting but low-dimensional example. Renormalization group studies of random 
magnets of Ref. [30] and techniques described in the review by ShankarEiJ must be considered to eventually understand 
the spin glass gap in presence of long-range Coulomb interaction too. 



C. Fermion propagator and spectral representation 



Perturbation expansions, which use the present model as a free (solvable) limit, are an interesting possibility to 
study for example itinerant spin glass models. In such expansions the fermion propagator, which is a site- localized 
propagator in time, should be known analytically. In particular its analytical properties are required in order to be 
able to evaluate diagrams of the perturbation theory. While it is not yet clear whether an exact analytical solution 
can be found in the ooi?5-B-limit, the present analysis allows to approach this solution in a qualitative way. The 
numerical study helps to find an analytical fit function, which could be considered as a diagrammatic element in the 
expansions mentioned above. Using our previous result for the density of states p(e) = —■^Im(G R (e)) the spectral 
representation 

/oo 
dup{u)/{iz n - u) (38) 
-oo 

with z n = (2n + l)7rfcsT ' jfi allows to evaluate the real part and thus the full Green function too. 
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FIG. 20. The real part of the fermion Green function is shown for 0- and 4RSB over a wide energy range. The zoomed 
region (upper left inset) shows all calculated 7?SB-orders (0, 1, 2, 3, 4). The result from Eq.^ using the OO-RS'-B-DOS model is 
contained (dashed curve carrying label e) and visibly different from 4RSB only in the zoomed inset. Gap edge singularities 
(spikes) with positions depending on RSB-order k are observed; curves are labeled by the order k. The logarithmic gap edge 
singularities are cut off in numerical calculations by small imaginary parts 8 = O(10 -6 ) in E + iS which ensure a finite distance 
from the branch cut of G. The free l/_B-behaviour is shown for comparison. 



Figure (inset) shows the step-by-step approach of the exact solution. The limit of spin glass interaction effects 
is shown by the crossover into the free limit l/E; together with this one observes an energy range where even the 
replica symmetric approximation is almost exact (roughly for > 3). The inset shows the increase in length of 
the well-approximated energy regimes with the growing order of RSB. The Figure also confirms the existence of 
hierarchical excitation energy regimes, where a certain minimal order of RSB guarantees an almost exact result. 
Moreover these regimes appear to scale with corresponding characteristic temperature regimes. We further observe 
that the logarithmic divergence of Re(G) at the Fermi level appears to be compatible for the scale chosen for the 
main Figure; in the zoomed inset the slope at Ep looks finite even in the interpolation model. This model however 
includes the logarithmic divergence which emerges only in an exponentially small (invisible) energy range near the 
Fermi level. 

The analytical result in QRSB for the retarded fermion Green function for T = and half-filling reads 



8 - A=± 



A exp{--c\ E )T(Q,-c\ iE ) - iirexp{-c\ E ){\ - erf(-i c x , E )) 



(39) 



where c\.e = E + Xx + i0+- The T-function in this expression reproduces the logarithmic singularities at the gap edge 
energies as seen in the numerical calculation above. Eq.^ crosses over into the free 1/iS-behaviour for large energies 
and is almost exact beyond \E\ > 3, ie for excitation energies more than four times beyond the spin glass temperature 
(recalling J = 1). From Eq.|38| it is clear that whenever p(u) drops to zero at a gap edge energy, the Cauchy principal 
value integral must diverge logarithmically. Only in the ooRSB limit is p no more discontinuous; it is even in u and 
no divergence occurs in Re(G(Q)). For finite Hubbard interaction the gap is not closed, but the density of states 
decays at least with some power at the Hubbard gap edge, which again prevents a singularity in Re{G). By means 
of the rapid convergence of RSB, the numerical analysis supports that the slope of G(E) at E — diverges. Only in 
case when the spectral density vanishes faster than linearly at E — the slope would remain finite. 



Figure 21 (fig21.gif) appended at e.o.f. 
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FIG. 21. Part a) shows Re(G(e)), e real, in 4RSB-approximation and superimposed the result obtained from the interpolation 
model for the exact density of states. Very good agreement is obtained. Both curves can hardly be distinguished except in the 
tiny regime near the small 4RSB-gap edge at |e| ~ 0.06, which is absent in the exact solution. Innocuous too, although present 
in the interpolating model solution in an exponentially small region near the Fermi level, is the logarithmically divergent slope 
of the real part of G(e). 

Part b) shows the effect of the Hubbard interaction for U — 1 (in units of J) in comparison with the U = 0-result contained in 
part a). For U = 1 the slope of Re(G) does no more diverge at the Fermi level. The Hubbard interaction squeezes a hard gap 
of size U into the spin glass pseudogap. The logarithmic divergence at the Fermi level is therefore removed for all finite U. 

In Fig.|l] part b) a solution for finite Hubbard interaction U = 1 is compared with the U — case; we justify the 
use of the ooi?5.B-model by the very good agreement with the calculated 4i?S*i?-curve (once the RS B-a,rtif&ct at the 
gap edge is omitted). The density of states together with Fig.^T] give a complete picture of the true (local) fermion 
propagator. Starting from the numerical data for arbitrary values of the Hubbard coupling (in units of the spin glass 
coupling J) we are ready to create, on the basis of computer algebra, objects that can be dealt with like standard 
Green functions and hence be used in diagram theories. 

The numerical analysis is completed by Figure ^2] for the purely imaginary Green function G(iz n ). Figure |22| shows 
a spectacular convergence already in the available low orders of RSB. Specific singular features, observed in G R (E) 
for real energies E, are absent on the imaginary axis (of course they are hidden and almost innocuous in G(iz n )) 
The absence of singularities on the imaginary energy axis gives rise to the speculation that an unexpected simple 
approximate form of the exact spin glass propagator may be good enough for calculations on the imaginary axis 
(leaving aside problems with the final analytic continuation to real energies). 
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FIG. 22. The (imaginary) fermion propagator G(ie) for half filling is displayed as a function of the fermionic frequencies 
z n = (2n+ Vj-RksT '/% — »t->o e which become continuous in the zero temperature limit. For (7 = 0, the dashed curve represents 
the 4_RS'_B-calculated one, which almost perfectly coincides with the one calculated from the DOS-interpolation model for full 
RSB. The Hubbard interaction effect is shown for the chosen strength U/J — 1. 



D. Scaling at the half-filling T = 0-transition 



The half-filling transition is sharp at T = and becomes continuous as x vanishes in the ooi?5B-limit. It can be 
driven by the variation of either e = E — [i = E — Ef which controls the pseudogap-shape, or by U — Ef which breaks 
down the central band. Finite temperatures act like a symmetry breaking perturbation which smears the transition. 
Our detailed i?5*i?-calculations suggest the following scaling behaviour 

p(e,T= 0,/x < U) « oi|e| + A a \e\ a gap regime (40) 

dC 

-f- ~ Inlel for e -> (41) 
de 

v - 1 - (fj, - U) 2 - 1 - q aa = 1 - qi for n > U (42) 
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where Fig.[lJ suggests a subleading exponent a — 2. Assuming that e scales like fj, — U near the transition one can 
propose the scaling form for the density of states 

p(e,n-U,T)=A\e\*f p (V—^,?-), (43) 

where, according to Eq.(40) our prediction is x = 1 for the infinite-range model, and f p is meant to be the scaling 
function of the density of states. 

For any finite fc-step RSB the transition is discontinuous; viewing the influence of smaller and smaller order parameters 
<7fc as k — > oo as a disopder fluctuation effect reminds of conclusions stating that disorder fluctuations can render a 
transition continuousc3'E3. The difference would be that we are concerned here with a one-dimensional effect in replica 
space, while those authors referred to real space fluctuations. 



IX. A SCENARIO FOR THE SPIN GLASS DRIVEN METAL-INSULATOR TRANSITION 



1. Non-half -filled system 



According to the preceding chapters, a central band always exists if the system is not half-filled. Then, the Fermi 
level must lie inside the upper spin glass pseudogap, which is separated from the Hubbard gap by the central band. 
A metal insulator transition must take place as fermion hopping increases beyond a critical value. Two possibilities 
arise: first the pseudogap defined by p(Ep) = survives until the metal insulator transition (MIT) takes place or, 
more likely, the pseudogap is gradually filled but states remain localized until the density of states becomes large 
enough. Then the MIT would take place. Its character must be different from standard Mott-Hubbard-Heisenberg 
metal insulator transitions, since the gap's existence rests exclusively on the frustrated magnetic interaction and, in 
addition, the localization of states stems again only from the randomness of this interaction. No further random 
potentials or random scattering is necessary to localize states in a way similar to Anderson localization. We explained 
already in the introduction how the random magnetic order might act in this way. In this context one should recall 
that broken time reversal allows for the unitary type of Anderson localization. In this respect it is an interesting 
detail that spin glass order is a special breaker of time-reversal invariance TRI, since the random magnetic moments 
break TRI locally, but globally the magnetization is zero without field and a global average breaking is not present. 



2. Half-filling 

In the half-filled system all gaps combine into one. At T = the spin glass gaps appear to be attached on each 
side of the Hubbard gap. The Fermi level lies in the center and one could expect the metal-insulator transition to 
be similar to that of the Hubbard model. However a detailed analysis must scrutinize this, in particular when the 
variance of the frustrated interaction is larger than the Hubbard coupling U. 
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FIG. 23. Expected weak filling of the spin glass pseudogap by fermion hopping in non half-filled systems; the shown form 
is assumed and not yet based on calculation. Anderson-like localization, generated by the random magnetic interaction, is 
indicated in the small density of states regime. The dashed line sketches the possibility of critical depletion at the metal-insulator 
transition as in the Coulomb case. 
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X. FERMIONIC SK-MODELS WITH HUBBARD-[/ DISTRIBUTIONS 



This paper has focused so far on effects of fixed Hubbard couplings and the final section can neither give a com- 
prehensive overview nor repeat with the same depth RSB- and Green function-calculation for [/-distributions. We 
limit the detailed analysis of random-[/ models to effects concerning the freezing temperature, including modifica- 
tions of tricritical behaviour We discuss related physical problems which necessitate further studies in the future. 
A lot of interesting possibilities emerge in models with random distributions of different Hubbard-couplings or with 
randomly placed Hubbard centers. Examples include long-standing physical problems like the ones mentioned in the 
introduction and described in Refs. IHH as well as new problems which appear in artificially tailored systems, or 
in mathematical toy models which help to probe the spin glass state by means of different kinds of [/-distributions. 
Apart from randomness, particularly interesting generalizations are sublattice structures modeling A-B alloys, for 
example with alternating signs of U, which allow staggered magnetic order in the sense of staggered deviations of the 
spin glass order parameter from its mean value. Such insulating states are nevertheless highly pregnant with compet- 
ing or perhaps coexisting magnetic superconductivity (or superconducting magnetism) if only hopping processes are 
taken into account. Ioffe and LarkinEj described disordered superconductors with smeared transition temperatures 
and randomly distributed couplings. They considered percolation of clusters, in which superconducting order can 
resist the destructive proximity effect of surrounding non-superconducting material. Below we consider a few differ- 
ent model classes which, according to what was achieved in the previous sections, shed some light on the interplay 
between randomness in U and the freezing temperature. Distributions of Hubbard couplings of the site-uncorrelated 
Wilson-Ising potential 

P({Ui}) =Af- 1 exp(-a 2 Uf - ajjf) 

facilitates the crossover from standard Gaussian distributions to a two level ^-distribution (random AB alloy). One 
may include site-site correlations by adding an additional Wilson-like term 

where either V denotes the discretized gradient on the lattice or one may think of the continuum limit. We will 
restrict the discussion to cases of a correlation-free distribution, which will not invalidate the chosen form of the spin 
glass order parameter. 



A. White noise distributed Hubbard couplings Ui 



We study first the case of a local Gaussian disorder correlation (white noise) described by 

((U i -(U i ))(U j -(U j )))=SU%j 
The [/-average of the replicated partition function considered by Eq.([|) can be written as 



(44) 
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(45) 



where TZ stands for the [/-independent part in the partition function of Eq.^| n and a are shorthand notations for 
charge and spin operators expressed in terms of Grassmann fields, which are integrated over by means of J 2?\& 
(expression n 2 — a 2 coincides with the square bracket of Eq.5). One can see that for Gaussian distributions what one 
needs to do is to replace U by U + SUxi and to perform the x^-averages at each site. 
The averaged freezing temperature can then be found from 
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(46) 



By the help of [/-averaged selfconsistent equations the following problems can thus be answered immediately: 

i) change in freezing temperature, 

ii) can continuous transitions occur down to zero temperature, 

iii) how strong is the depression of spin glass order and of Tf by attractive interaction (negative U) in comparison 



26 



with the restoring repulsive interaction (positive U); distributions can be defined which probe this competition by the 
relative weight, the limit of infinitely broad Gaussian {/-distributions being a special case, which should render the 
mean coupling (U) irrelevant, and 

v) can random alloy models with special permitted (C/)-values be expected to exhibit significant behaviour? 

According to Eq.[h| SU cannot be absorbed in a /i-shift. The asymmetric [/-dependence (attractive U suppress 
spin glass order while repulsive interactions support it) allows a crossing of the critical curves Tf(5U). As one can 
further observe in Fig.|24|, the critical temperature Tf (SU) averages over the competing effects of positive and negative 
U and varies little for broad distributions. 

U=l 




8U 2 5 

FIG. 24. The variation of the freezing temperature Tf under increased variance SU of the Gaussian [/-distribution is shown 
for three chosen values of the Hubbard coupling, U = — 1 (attractive), U — 0, and U = 1 (repulsive). For each U, Tf(5U) 
is shown for three different chemical potentials (fi- value given at the right end of T/-curve), which are symmetrically placed 
around the half-filling condition n = U/2 for SU = 0. 

The selfconsistency equation for Tf is simplified for SU — > oo and reads 



Tf = l/ 



2 + e (- M -l/2)/7> 



(47) 



which is U -independent because of the infinitely broad distribution. The freezing temperature approaches Tjjjx) = h 



for /i — > oo but drops discontinuously to zero near the stability limit at /i s ^ 

. 5 F 



-0.7315 as shown in Figure 
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FIG. 25. Solution of 7/ in the limit of infinitely broad Gaussian distribution of U : one finds Tf = \ in the limit of large 
fi — * oo while, in the small /i-regime, a stability limit shows up at fj, s .i. — —0.7315 and implies a discontinuous drop to Tf = 0. 



The limits SU — > oo and /i — > oo do not commute. For fixed [/-variance and arbitrary U one finds that the freezing 
temperature decreases exponentially for /i — > oo. 
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For several differently chosen distributions the averaged Ty-equation results in Figure 26 showing an exponential 
decay in the large /x-limit of the critical temperature (for fixed U). A universal temperature T/o = 7/(mo) is observed, 
where Tjq is independent of distributions which are symmetric with respect to U. The value /io of this invariant point 
depends linearly on U and obeys Ho(U) — 0.9939(L/ + 0.8626). One may check that the critical temperature Tf given 
by Eq. ^ and Figure [l] at agrees too. 




FIG. 26. Exponential decay of Tf for the normalized ^^distributions a) ~ exp(—x 2 /2)/(x 2 + 1), b) ~ l/( a;2 + 1), c ) 
■ exp(-x 2 /2), d)- exp(x 2 - x i /2), and e) ~ exp(x 2 - x 4 /10) with fixed U(= 1). 



For negative /i the freezing temperature vanishes discontinuously. A line of tricritical points exists. The exponential 
tails at large positive \i originate in the rare large [/-regions. None of the given examples allows a continuous phase 
transition at T = 0, since Tf remains nonzero for all finite fi (allowing only a discontinuous drop to TV = 0). 



B. Alloy models with two different Hubbard couplings 



Random alloys with only two different Hubbard couplings Ui, U2, which are assumed to be realized with probability 
w\ and W2 = I — Wi respectively on each site obey 

P{U t ) = Wl S(U t - U x ) + (1 - w 1 )S{U i - U 2 ) (48) 

In this model the spatial homogeneity of the spin glass order parameter (Q) remains unaffected. We observe that the 

Tf-solution for this model also passes through the universal point of Figure |26| provided W\ = | . 

A highly interesting AB-a\\oy model with staggered Hubbard-interaction, defined with a [/-distribution 

P({Ui}) = S(Ui - U A )5 l , lA + 5(Ui - U B )8 itiB , (49) 

can break the spatial homogeneity of (Q). In case that V a on sublattice A is strongly negative, this sublattice will 
become nonmagnetic, while the range-free magnetic interaction maintains glassy order on sublattice B provided Ub 
is repulsive for example. When nearest neighbor hopping of sufficient strength is introduced, an extreme proximity 
effect between superconductivity and magnetism will result: local pairs can only become delocalized when single 
fermions tunnel through the magnetic sites. A strong filling dependence must be expected. In addition the strong 
interference of an eventually superconducting sublattice system with glassy magnetic order is a challenge for research 
on microscopic superconducting glass phases. 



XI. ANTIFERROMAGNETIC- AND FERRIMAGNETIC GLASSY PHASES 



The simplest model with tange-free magnetic interaction, which nevertheless supports antiferromagnetic order, is 
the Korenblit-Shender modeEJ. The magnetic interaction is supposed to act only between A and B atoms, which one 
may imagine to sitjan neighbouring sites, with a mean antiferromagnetic form together with random fluctuations. 

In previous workEa we solved this model in one-step RSB and generalized it in several directions. Spin glass order 
parameters Qa, Qb within each sub-lattice differ in the glassy antiferromagnetic phase transition, and a third field is 
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necessary to describe the A — B coupling. In a field we found a level crossing of elements of Qa and Qb, which led us 
to define a new i?S*i?-scheme. This scheme is allowed to be A — B asymmetric which, in the strongly A~ B- asymmetric 
ferrimagnetic phase, turns out to be better than the standard solution. Details will be published elsewhere. 
The Hubbard interaction can be built into the antiferromagnetic model in the same way as described before. Hubbard 
interactions Ua, U b of different strength can enhance the A-B asymmetry. This is evident from the fact that sufficiently 
negative U, say on A-sites, render these sites nonmagnetic, while sites with positive U, say on i?-sites, maintain glassy 
order. This case has very strong asymmetry and the new RS B-scheme of Ref. |l8| must be considered as a candidate 
to describe such phases. 

Such problems are left for future research. When hopping is also allowed for, superconductivity on the A-sublattice 
in competition with glassy magnetism on £?-sites revive the problem of microscopic superconducting glasses as well. 

XII. SUMMARY OF RESULTS AND OUTLOOK 

We presented the analytical solution for the generating functional of a model with competing Hubbard- and frus- 
trated range-free Ising spin interaction together with detailed numerical evaluations of the resulting self-consistency 
equations in order to describe relations between magnetism and electronic properties. The model was interpreted 
as the localized (but non-local) limit of a Hubbard model with additional frustrated spin interaction. The discrete 
time slicing method was, to our knowledge, used for the first time to solve a spin glass problem analytically (all inte- 
grations were performed before the continuum time limit was taken). Taking advantage of the rapid convergence of 
/c-step replica symmetry breaking and using exact relations we drew conclusions on the ooRSB solutions of important 
quantities such as the fermion propagator G(E). The 2 — 4-step approximations were found to be close enough to the 
exact solution (by means of comparison with exact relations) and allowed us to identify and to eliminate artifacts of 
finite-step RSB in the extreme low temperature regime. By comparison of 0-,l-,2-,3, and 4-step results, we identified 
hierarchical energy regimes, < |e| < §P°\ where corrections beyond a certain fc-step-approximation are negligibly 
small. For decreasing excitation energies above or below the Fermi level an increasing number of steps is required at 
T = for good accuracy. 

We found that all T = phase transitions of the model (including its generalization to random distributions of 
the Hubbard coupling) are discontinuous with the exception of half- filling transitions. Those become continuous in 
the limit of infinite-step RSB (i.e. their discontinuities decay with increasing number of RS £?-steps and vanish in the 
limit ofooRSB). 

The density of states around the spin glass pseudogap at Ep was found to obey the linear dependence p(E) = 
c\\E— Ep\ with ci ss 0.3 in a wide gap regime and c\ 0.13 asymptotically close to the Fermi level. This was concluded 
from a fast convergence of p'(E) under increasing number of RSB-steps. The linear behaviour of p reminded us of the 
Efros-Shklovskii Coulomb-gap in two dimensions. These gaps have totally different origin: the Efros-Shklovskii gap 
originates in the bare Coulomb interaction (and naturally involves a dependence on space dimension) while the present 
one is caused by the spin glass order. In finite-range models the latter is expected to depend on space dimension 
predominantly through the destruction of glassy order at its lower critical dimensions for finite temperature and 
T = 0-transitions respectively. 

We believe that a renormalization group analysis acting in replica space (not in real or momentum space) will be 
helpful to describe the asymptotic behaviour linked to the hierarchical energy scales, which were discussed in our 
paper for the range-free, model. Eventually a new time of RG must be created for this purpose. In accordance with 
general considerations!!!!, Johnson and KhmelnitskiO showed for a disordered system with long-range (unscreened) 
Coulomb interaction how important and useful a RG-technique can be for the understanding of pseudogaps. 

From our numerical analysis we predicted a set of other scaling laws. The linear decay of the density of states led us 
to claim a divergent slope of the Green function's reaL-oarf at Ep and hence a divergence of a corresponding vertex, 
which results from charge conservation Ward identityH. These divergences are logarithmic and become visible only 
exponentially close to the Fermi level. This reminds of the extremely slow approach to equilibrium, which is usually 
discussed for spin dynamics. In the present case, in the absence of spin dynamics it is the quantum fermion dynamics, 
generated by the noncommuting fermion- and Hamilton operator, which experiences the hierarchical i?5-B-structure of 
spin glass order parameters. For the fermion concentration (integrated density of states at T = 0) and for the related 
T —>■ limit of the integrated Parisi function, qi — lim^^o f dx q(x,T), we predicted v—1 ~ 1 — q\ ~ (p— U) 2 , p > U , 
near the half-filling transition. The transition generates an important crossover of the Fermi level from the center 
of a spin glass pseudogap into the Hubbard gap and hence into a deeper insulating state. We stressed its role for 
derealization in itinerant models and described related scenarios for metal insulator transitions. 

The present approach was intended to progress towards an analytical solution at T = and the construction o£-a 
low temperature expansion. A generalization to fermionic spin glasses of a recent numerical study by Crisanti et aO, 
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derived from an optimization algorithm and analogy with T — standard spin glasses, could be a new alternative 
route way as well as the application of a renormalization group in replica space. 

The models we studied discarded the long-range part of the Coulomb interaction. This interaction, which is known 
to produce the Efros-Shklovskii Coulomb pseudogap, should be taken into account in a more complicated model. The 
eventual combination with the spin glass pseudogap should be analyzed. 



We finally note that the smooth form obtained for the fermion propagator in section |VIII Cj , and its analytical 
approximation, can be a starting point of a perturbation theory) or at least control limit in studies of more complicated 
cases such as itinerant models with fermion hopping, mixed valence (Anderson-)models, and s-d coupling models 
between spins of mobile carriers and spin glass ordered localized spins. 
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